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This set of exercises relates to Book I of SM358 and can be used for further practice and for revision purposes. To 
gain maximum benefit, it is important to attempt each question for yourself before looking up the solution. 


Exercises 
Topic 1 — de Broglie waves 
Exercise 1.1 A one-dimensional de Broglie wave describing a free particle takes the form 
Wap (x, t) = eve). 
where A = 3.20 x 10!9m~!, B = 7.44 x 10!4 s7! and the potential energy is equal to zero everywhere. 


(a) What are the wave number, angular frequency, wavelength and period of this wave? 


(b) What are the momentum magnitude, kinetic energy and mass of the particle? 
Exercise 1.2 Show that the one-dimensional de Broglie wave function 
Wap (x, t) = eu). 


where A and T are positive constants, is an eigenfunction of the momentum operator P}, and find the 
corresponding eigenvalue. Is your answer consistent with the de Broglie relation linking the magnitude of 
momentum to the wavelength of a free particle? 


Topic 2 — Operators 

Exercise 2.1 Determine whether each of the following operators is linear: 
O1f (2) = 2? f(z) 
Oof(x) = f?(x) 
Osf(x) = f(z + a), 


where a is a constant, and f(x) is an arbitrary function. Explain your reasoning. 
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Exercise 2.2 Simplify the expression 
(0102 = O201)g(2) 
where the operators Ô and O» are defined by 


ôs) = YO and Oof(x) = x° f(x) 


for any function f(x). 


Topic 3 — Schrödinger’s equation 


Exercise 3.1 An isolated one-dimensional system consists of two particles of masses mı and mo interacting 
via a potential energy function 5C (x1 — x2)”, where C is a positive constant and x, and x2 are coordinates 
representing the positions of the particles. 


Write down a Hamiltonian function, a Hamiltonian operator and Schrodinger’s equation for this system, showing 
explicitly all the variables and parameters involved. 


Exercise 3.2 A one-dimensional system consists of two particles of masses mı and m2, which do not interact 
with each other, and are each subject to a potential energy function of the form 5Cx?, where C is a positive 
constant and x represents position. 


Write down a Hamiltonian function, a Hamiltonian operator and Schrodinger’s equation for this system, showing 
explicitly all the variables and parameters involved. 


Exercise 3.3 (a) Write down Schrodinger’s equation for a free particle of mass m in one dimension and show 
that the solutions of this equation obey the principle of superposition. 


(b) A beam of electrons is incident on an absorbing screen which contains two narrow slits. Some particles in the 
beam pass through the slits and are detected on a screen on the far side. How does the principle of superposition 
help to explain the pattern built up by a large number of particles falling on the detecting screen? 


Topic 4 — Separation of variables 


Exercise 4.1 The one-dimensional form of Schrodinger’s equation for a particle of mass m in a potential energy 
well V(x) is 


h? o? 


a, 


T W(x, t) + V(x) WV (2, t). 


(a) Use the method of separation of variables to show that solutions to this equation exist in the form 
(x,t) = w(x)T(t). Derive the equations that y(x) and T(t) must satisfy, and solve the equation for T(t). 


(b) Can all the solutions of Schrodinger’s equation be written in this form? Explain your answer. 


Exercise 4.2 The time-independent Schrodinger equation for a particle of mass m inside an infinite 
three-dimensional square well takes the form 

h2 82 82 32 

-— | =— t = +t = =E ; 

2m = Tr Oy? 4 Z) W(x, y, z) p(z, y, z) 
Use the method of separation of variables to show that solutions to this equation exist in the form 
y(x, y, z) = X(x)Y (y)Z(z). Derive the equations that X(x), Y (y) and Z(z) must satisfy, and specify a 
relationship that exists between the constants appearing in these equations. 


Topic 5 — Born’s rule 


Exercise 5.1 (a) State Born’s rule for a wave function Ẹ¥ (x,t) of a particle in one dimension. 


(b) A normalized one-dimensional wave function takes the form 


[9 
Viet) = T sin (=) e HZt/h foro <t ÁL 


where L is a length and F is an energy. Estimate the probability that a measurement of the particle’s position will 
give a value in the small interval between x = 0.249L and x = 0.251 L. Does this probability depend on time? 


Exercise 5.2 At time t = 0, a particle is in a state described by the normalized wave function 


2 a: 51x 

£ sin (7 forO<a< L, 
U(x,0) = L ( L ) eS 

0 elsewhere. 


What is the probability of finding the particle somewhere between x = 0 and x = L/5 at t = 0? 


You may use the standard integral 
© a T 
du = >. 
f sin udu = > 


Exercise 5.3 A particle is in the ground state of a one-dimensional finite potential well with walls at x = — L /2 
and x = +L /2. Outside the well, the ground-state energy eigenfunction is 


Ce" forg > +L /2, 
yı (x) E +ax / 
Ce for x < —L/2. 


(a) What is the probability of finding the particle somewhere outside the well? Express your answer in terms of 
C, a and L. 
(b) What is the probability of finding the particle somewhere inside the well? 


Exercise 5.4 At t = 0, a particle is described by a one-dimensional wave function 
Y(x,0) = Az en 
where A and a are constants. Where is the particle most likely to be found? 
Exercise 5.5 (a) Use Born’s rule to show that the SI units of a wave function in one dimension are in 


(b) Ifq@isa length and z is a position coordinate, which of the following functions of x are unsuitable, on 
dimensional grounds, for representing a one-dimensional wave function at any given instant in time? 


(i) f(z) = e772? /2a? 


(ii) g(x) — Ta 
a 

is _ m —z?/2 

Gii) A(z) =4/ a : 


Topic 6 — Normalization and phase 


Exercise 6.1 Att = 0, astationary-state wave function takes the form 
2/92 
U(r,0) = Aew® 2”, 


where A is a constant and a is a positive constant. Use the normalization condition to find a possible value for A, 
expressed in terms of a. How much freedom is there in choosing this value? 


You may use the standard integral 


OO 
f eo /@ du=a/n fora> 0. 
—oo 


Exercise 6.2 An unnormalized one-dimensional energy eigenfunction takes the form 


sin (=) forO<a< JL 


0 elsewhere. 


p(z) = 


Obtain the corresponding normalized eigenfunction (defined up to a phase), explaining your reasoning. 


You may use the standard integral 


nT 
è NT P 
f sin? u du = Ea for integer n. 
0 


Exercise 6.3 An energy eigenfunction for a particle in a three-dimensional infinite square well takes the form 


W(x, y, z) = Asin (=) sin =) sin (=) 


forO<a2<L,0<y< Land0 < z < L, with y(x, y, z) = 0 elsewhere. Find a suitable value of A that ensures 
that this eigenfunction is normalized. Explain your answer. 


You may use the standard integral 


nT 
x NT : 
/ sin? udu = a for integer n. 
0 


Exercise 6.4 In this question, w, and Yz are two different normalized energy eigenfunctions and 


_ pi + ie 


Which of the following combinations of Yı and Yə (where a is real) represent exactly the same state as Ya? 


VA 


_ 1 — ip 
vp = -A 
ip -y2 
Ye = E 
_ ipi- yp 
vp = = 
_ ipi ty 
pe = -A 
Je- eit + eila tT/2) yo 
o v2 
bo = ella—m/2) qh, iil el aby 


V2 


Topic 7 — Square wells 


Exercise 7.1 A particle of mass m is in the ground state of a three-dimensional infinite square well which 
occupies a cubic region with sides of length L. By how much does the energy of this system change if the sides of 
the cubic box are halved in length to L/2, with the particle remaining in the ground state of the new box? 


Exercise 7.2 A spinless particle of mass m is in a three-dimensional infinite square well which occupies 
a cubic region with sides of length L. How many degenerate quantum states correspond to the energy level 
E> = 3h? /mL?? 


Topic 8 — Continuity boundary conditions 


Exercise 8.1 Write down continuity boundary conditions for an energy eigenfunction w(x) and its derivative 
dy)/dx at the walls of a one-dimensional finite square well. 


In one or two sentences state how the behaviour of an energy eigenfunction at the walls of a finite square well 
differs from the corresponding behaviour at the walls of an infinite square well. 


Exercise 8.2 A particle is in the ground state of a one-dimensional finite square well with walls at xv = —L/2 
and z = +L /2. Inside the well the energy eigenfunction describing the particle is 
Yin(z) = Acos(kzx), 
where A and k are constants. Outside the well, for x > +L/2, the energy eigenfunction is 
Pout) =Ce™, 
where C and a are constants. Apply continuity boundary conditions at x = +L/2 to find a relationship between 


k, L and a that does not include the constants A and C. 


Exercise 8.3 A one-dimensional infinite well is described by the potential energy function 
0 for—L/2<2<0 
VYaie Vo for0<2<L/2 
co elsewhere, 
where Vo is a positive constant. An energy eigenfunction for a particle inside this well takes the form 
Acos(kiz)+ Bsin(kiz) for —L/2<x<0 
p(z) = C cos(kax) + Dsin(kgx) for0 < x< L/2 
0 elsewhere, 


where A, B, C, D, kı and kə are constants. 


Use the time-independent Schrödinger equation to find a relationship between k1, k2, Vo and the mass m of the 
particle. 


Exercise 8.4 Consider again the situation described in Exercise 8.3. Use appropriate boundary conditions at 
x = —L/2, x = L/2 and x = 0 to find a relationship between k1, k2 and L that does not involve the constants A, 
B,C or D 


Topic 9 — The overlap rule 


Exercise 9.1 (a) State the overlap rule for energy measurements in a state described by the wave function 
(x,t) in a system with discrete non-degenerate energy levels. 


(b) Supposing that a system is described by a wave function that is an even function of x, use the overlap rule to 
show that there is no chance that any measurement of the system’s energy will give a value that corresponds to an 
energy eigenfunction that is an odd function of x. 


Exercise 9.2 A system has energy eigenfunctions w(x), W2(x), ..., with corresponding energy eigenvalues 
E1, E2, .... If the system is in a state described by the energy eigenfunction w(x), what is the probability of 
getting the value Hy Æ Ei in an energy measurement? Use your answer to show that the overlap rule implies that 
the eigenfunctions ~1(x) and ~2(x) must be orthogonal. 


Exercise 9.3 A particle of mass m is in a one-dimensional infinite square well with walls at x = —L/2 and 
x = +L/2. At t = 0, the particle is described by the normalized wave function 


W(e,0) = z for —w/2 < x < w/2 
0 elsewhere, 


where 0 < w < L. If the energy of the particle is measured at time t = 0, what is the probability that the 
ground-state energy eigenvalue E = 72h? /2mL? will be obtained? 


Topic 10 — The coefficient rule 


Exercise 10.1 At time t = 0, the state of a particle is represented by a wave function of the form 


W(x, 0) = Zvi (2) + ayola), 


where pı (a) and w2(x) are normalized energy eigenfunctions corresponding to different energy eigenvalues Fy 
and Ev. Find a possible value of a that ensures that U(x, (0) is correctly normalized. Would the state be changed 
by changing the phase of a? 


Exercise 10.2 At time t = 0, the wave function describing the state of a system takes the form 
W(x, 0) = 2Ny1(x) — V5Ny0(z), 


where w(x) and Y2(x) are normalized energy eigenfunctions with energy eigenvalues FE, = £ and Ey = 2e 
respectively. Find a suitable value for the normalization constant N, and determine the expectation value of the 
energy at time t = 0. 


Exercise 10.3 A system is in a state described by the wave function 
W(x, t) = A(dy (we Ft — ino (aye ith), 
where w(x) and Y2(x) are orthonormal energy eigenfunctions corresponding to the eigenvalues Æ; and E 


respectively. Find a suitable value for the normalization constant A and determine the expectation value of the 
energy at any time t. 


Exercise 10.4 Three normalized stationary-state wave functions y1 (x) eE», p(x) e'2"/" and 


palz) e iBst/ i corresponding to energies /, E2 and E3, are added with coefficients a1, az and a3 to produce the 
wave function 


W(x, t) = a yi(x) ei + a a(x) eFt + ag da (x) eitr, 
(a) Write down the condition that must be satisfied by a1, a2 and ag in order for the wave function U(x, t) to be 
normalized. 


(b) If the coefficients a1, ag and ag are all multiplied by —1, will this make any difference to the expectation 
value of the position x? 


(c) If a, alone is multiplied by —1, will this, in general, make any difference to the expectation value of the 
position x? 


(d) Ifa alone is multiplied by —1, will this, in general, make any difference to the expectation value of the 
energy, E? 


Topic 11 — The sandwich integral rule 


Exercise 11.1 (a) State the sandwich integral rule for the expectation value of an observable A in a system of 
one particle in one dimension. 
(b) At t = 0, a particle is in a state described by the wave function 

W(x,0) = Ae, 


where a and A are constants and a > 0. Find the expectation values of x and py in this state. 


Exercise 11.2 The normalized energy eigenfunctions of a particle in a one-dimensional infinite square well with 
walls at x = 0 and x = L are 


abn (2) f= sin (=) forO<a< L 
nv) = 


0 elsewhere, 
form 1; 2. ssa 
Calculate (x°) in the state with quantum number n. 


You may use the standard integral 


nT 
I u? sin? udu = Mt (2n?r? — 3). 
A 12 


Exercise 11.3 Find the expectation value of the observable p2 in the state described by the normalized wave 
function 


1 1/2 _2? /2a? 
wiz, 0) = We e ; 


You may use the fact that 


d? 2/972 1 (x? 275,2 
Baa —a*/2a ene (ee —a*/2a 
dx? (e ) a? (5 1) ° 


together with the standard integrals 


f e” du = yr and | we du = A 


T 2 


2y 


Exercise 11.4 Find the expectation value of the observable e~* ? where À is a real constant in the state 


described by the normalized wave function 


1 ee 2? /20? 
Y(zx,0) = Sie e : 


You may use the standard integral 


co 
oi T 
f e * du = 4/—, 
_ Va 


for any positive constant, a. 


Topic 12 — Uncertainties and the uncertainty principle 


Exercise 12.1 At t = 0, a particle is in a state described by the normalized wave function 


¥(2,0) = yila) + iyala), 


where 1 (x) is a normalized energy eigenfunction with eigenvalue Ey = £ and 72(2) is a normalized energy 
eigenfunction with eigenvalue E = 2e. Find the expectation value and uncertainty of the energy in this state. 


Exercise 12.2 A particle is in a state characterized by the expectation values (x) = 1.0nm and (x?) = 2.0 nm?. 
(a) What is the uncertainty in x in this state? 


(b) According to the Heisenberg uncertainty principle, what are the minimum values of the uncertainties in pz 
and p, that could be consistent with this uncertainty in x? 


Topic 13 — Harmonic oscillators 


Exercise 13.1 Verify that the function ¢)(x) = e722", where a is a positive constant, is a solution of the 
time-independent Schrodinger equation for a harmonic oscillator: 


are + 5Cx*y(z) = Ey({x), 


provided that a has an appropriate value, which you are asked to find. Also find the corresponding energy 
eigenvalue, E. 


Exercise 13.2 A harmonic oscillator has a zero-point energy of 1.5eV. What is the wavelength of the light 
emitted when this oscillator makes a radiative transition from a state with quantum number n = 5 to a state of 
lower energy? 


Exercise 13.3 Use raising and lowering operators to show that 
RS ih 
in any energy eigenstate of a harmonic oscillator. 


Exercise 13.4 (a) If A and Al are lowering and raising operators for a harmonic oscillator, show that 
OO 
jz f va(a)[(A +i + (A- A") n(x) de =, 
—Co 


where w,,(a) is an energy eigenfunction of the harmonic oscillator. 


(b) Hence use the formulae 


t= (A+8') and P, = — 


to show that 
2 2 2 
(x) _ a (pz) 
ae h2 
in any energy eigenstate of the harmonic oscillator. Given that at = h?/mC, where m is the mass of the 
oscillating particle and Cis the force constant of the harmonic oscillator, what is the physical interpretation of this 
result? 


Topic 14 — Bound-state wave packets 


Exercise 14.1 At time t = 0, a wave packet describing the state of a particle in a harmonic potential energy well 
takes the form 
1 ; 
Y(x,0) = D (Yıl) + ip2(z)), 
where Yı (x) and Y2(x) are normalized energy eigenfunctions of the harmonic oscillator with eigenvalues 3 hwo 
and 3hwo, respectively, and wo is the classical angular frequency of the oscillator. 


(a) Write down an expression for U(x, t) at any later time t, expressing your answer in terms of wo. 


(b) Are there any times after t = 0 when Y (x,t) is orthogonal to the initial wave packet, U(x, 0)? 


Exercise 14.2 (a) A wave packet describing the state of a particle in a harmonic potential energy well takes the 
form 

Jy (Pole) eM" + ala) eH"), 

where wo(x) and Y2(x) are normalized energy eigenfunctions for the ground state and second-excited state of the 
oscillator, with energy eigenvalues Eo = hwo and E2 = 3hwo, respectively, where wo is the classical angular 
frequency. Show that the probability density for this wave packet oscillates with period 7 /wo. 


T(z,t)= 


(b) What is the relationship between Y (x,t) and ¥(—zx, t)? Does the motion of the wave packet correspond to a 
to-and-fro classical oscillation? 


Topic 15 — Free-particle wave packets 


Exercise 15.1 At time t = 0, a one-dimensional wave packet for a free particle of mass m takes the form 


1 i ; 
WU(x,0) = = A(k) ei** dk. 
(0,0) = = | AQ) 
Write down an expression for the wave packet at a later time, t, defining any quantities you introduce. What is the 
physical significance of the function A(k) with respect to measurements of momentum? 
Exercise 15.2 (a) At time t = 0, a one-dimensional wave packet for a free particle of mass m takes the form 
1 aiaz 
=e for —L/2 < x < L/2 
W(x,0) = 4 7 
0 elsewhere, 
where q is a real constant. What is the momentum amplitude function for this wave packet? 


(b) Assuming that a = 1/Z and using the fact that 


A estan 2 
f sin (4) ay = 0.897, 
0 u 


calculate the probability that a momentum measurement taken on the wave packet at t = 0 will give a value 
between —fa and +ħa. 


Exercise 15.3 (a) Explain why the following solution to the free-particle Schrödinger equation: 
Wap (x, t) = Acie 24/2), 


where A is a constant and Ej, = h?k?/2m, is not entirely satisfactory for representing the state of a free particle in 
one dimension. 


(b) Consider the following expression for a normalized free particle wave packet: 


1 ie 

—— | A(k) el? Fet/Mag, 

V 2r j ( ) 

Write an expression involving A(k) that is equal to the probability that a momentum measurement made on a 
particle described by U(x, t) will yield a value between ik, and hko. 


P(x, t) = 


What condition must A(k) satisfy in order for Yag (x, t) to be normalized? 


Topic 16 — Ehrenfest’s theorem 
Exercise 16.1 State Ehrenfest’s theorem for (x) and (py) for a particle of mass m subject to a potential energy 


function V(x). Show that, if V(x) = mga, where g is a constant, then Ehrenfest’s theorem gives 


d’) _ 
dt? 


Topic 17 — Scattering and tunnelling 


Exercise 17.1 A mono-energetic beam of particles, travelling in the x-direction is incident on a step located at 
x = 0. In Region 1 (x < 0) the beam is described by the energy eigenfunction 


i (a) _ Aei*1# she Be7ik2 
where A and B are constants. In Region 2 (x > 0) the beam is described by the energy eigenfunction 
polz) — Cel? 


where C is a constant. Use appropriate boundary conditions at x = 0 to relate A and B to C and hence determine 
the ratio B/A in terms of the wave numbers kı and kə. 


Exercise 17.2 A mono-energetic beam of particles, travelling in the x-direction is incident on a step located at 
x = 0. In Region 1 (x < 0) the beam is described by the energy eigenfunction 


wy (x) — Aei*1# + Beir, 
where A and B are constants and k; is real. In Region 2 (x > 0) the beam is described by the energy eigenfunction 
pale) = Ce, 


where C is a constant and a is real. Use appropriate boundary conditions at x = 0 to relate A and B to C and 
hence show that |B/A|? = 1. 


Exercise 17.3 Find the probability current associated with the free-particle energy eigenfunction y(x) = Ae‘**, 
where A is a constant. 
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Solutions 
Topic 1 — de Broglie waves 


Solution 1.1 (a) A one-dimensional de Broglie wave, travelling in the positive x-direction, takes the form 
Yag (x, t) E eee. 


where k is the wave number and w is the angular frequency of the wave. For the de Broglie wave given in the 
question we have 


b=A=320x10" m and w= B = 7.44 x 10% =: 

The wave number and angular frequency are related to the wavelength À and period T by k = 2r /AÀ and 
w = 2r /T, so 

27 27 10 
= — = ~ = 1.96 x 10 

k  3.20x 10m z 

and 

20 27 


= = -15 
w  7.44x10l4s-71 © 8.45 x 107" s. 


r= 


(b) The momentum magnitude is 

p = ħk = 1.06 x 107°4J s x 3.20 x 101? m™t = 3.39 x 107° kg ms™t. 
The energy of the particle is 

E = ħw = 1.06 x 1074 J s x 7.44 x 1014 s = 7.89 x 107” J, 


and, since the particle is free with zero potential energy everywhere, this is equal to the kinetic energy, Ekin- 


Finally, the kinetic energy and momentum magnitude of the particle are related by Ein = p?/2m, so the mass of 
the particle is 
p? (3.39 x 10724 kg m s71)? 


= = = 7.28 x 10° kg. 
= OE 2x 7.89 x 10-2J E 6 


Comments (1) The units of À, T and p are all correct. For mass we have J = kg m? s7?, so 


(kgms™!)?/J = kg, as required. 


(2) You may know from elsewhere in physics that the phase speed of any wave is given by Uwaye = fA = w/k. 
Please note, however, that the phase speed of a de Broglie wave is not equal to the speed of the free particle it 
describes, which is given by p/m. According to the last equation in the solution, the speed of the particle is 
p/m = 2Exin/p = 2w/k, which is twice the phase speed of the wave! 


Solution 1.2 We have 
~ ð ; 
Pr Vas (£, t) = i RA 


— ih 27) e2rile/à—t/T) 
AÀ 
_ 20h oni(a/d—t/T) 
X ? 


so Vap(z, t) is an eigenfunction of p,., with eigenvalue 27i/A = h/A, where h is Planck’s constant. The positive 
constant \ is the wavelength of the wave (i.e. the distance between points that have the same phase at a given 
instant). So this result is consistent with the de Broglie relation, \ = h/p. 


Comment To check whether any function is an eigenfunction of a given operator, apply the operator to the 
function and see whether you get the same function back, multiplied by a constant (the eigenvalue). 
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Topic 2 — Operators 
Solution 2.1 O; is linear because, if A and yu are arbitrary constants, 
O1 (AF (2) + n9(x)) = 2? (Af (a) + ug(2)) 
= Ax? f(x) + ux? g(x) 
= \(01f (x) + n(Org()). 
Ô» is not linear because, if A and p are arbitrary constants, 
Oo(Af(x) + ug(2)) = [AF (2) + ugla)? 
= f(a) + p’g’ (x) + 2Au f(x) g(2) 
= X? (Oof(x)) + 1? (Oog(a)) + 2du f(x) g(x) 
# (Oof(x)) + u(Ozg(x£)) Gn general). 
Ô; is linear because, if \ and y are arbitrary constants, 
Os(Af (x) + ug(2)) = Af (£ + a) + pg(a + a) 
= \(Osf(x)) + u(Osg(2)). 


Comment To check whether a given operator is linear, we always start by applying the operator to a general 
linear combination of functions Àf + jg and then try to rearrange the result. 


Solution 2.2 We have 
01 029(x) = O: (O29(z)) 


while 


dg(x) 
dx ` 


02019(z) = O2(O19(2)) a 
Subtracting these two results we obtain 


(0,02 = 0201) 9(z) = 2gg(x). 


Comments (1) When simplifying 0,0. (g(x)), note that the derivative operator 0; = d/dz acts on x2g(x), 
which must be differentiated using the product rule. This gives a different result from 0201 9(x) = x7dg/dz, so 
the order of the operators matters (they do not commute). 


(2) Since the final equation is valid for any function g(x), we can also express it as a relation between operators: 
0102 — 020; = 22, 


where the expression on the left-hand side is called the commutator of the operators O1 and Oo. 
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Topic 3 — Schrodinger’s equation 
Solution 3.1 The Hamiltonian function is 
2 
Pr 


2 
H= DE Pr,2 
2m1 2Moə 


+ $C(21 — 2)’, 


where pz, ı and pz 2 are the momenta of the two particles. 
Since P = —ihO/Ox and P? = —h?9?/Ax?, the corresponding Hamiltonian operator is 


~ ho o h? 6? 
head o a aaa 
2m1 Ox? 2m Ox > 2 (2i 2) 


Schrodinger’s equation then takes the form 


ð h? oy t h? OU («x t 
ihR—V (21, x,t) 2 (zL 22, ) oe (1,225 ) 

Ot 2m oxi 2mMə x5 
Comment We introduce separate position variables (xı and x2) and momentum variables (pz, and pz,2) for the 
two particles. The wave function Ẹ (x1, £2, t) of the two-particle system depends on both xı and z2, as well as on 
time, t. 


Solution 3.2 The Hamiltonian function is 


2 2 
Pot a Pr,2 


= 2m 2me 


+ 30r + 5Cx5. 


The corresponding Hamiltonian operator is 
re o io 


Heo —— hd + lOr. 
2m; ðx? moa Tit geta 


Schrödinger’s equation then takes the form 


o h? PU (x1 “2 t) hi? OU (x1 T2 t) 
ih—wU . t = e] , ’ e] 
i (1, 22,4) 2m4 dx? Ime ax? 


a + (40x? + 4Cx}) (x1, £2, t). 


Comment In Exercise 3.1 the particles interacted directly with one another via a potential energy function that 
depends on (a1 — x2)”, the square of the separation of the particles. The situation described in this question is 
different because the two particles do not interact with one another, but each particle is subject to a potential 
energy function (due to an external agency) which depends on the square of the separation of the particle from the 
origin (5Cxz for particle 1 and $C x3 for particle 2). 


Solution 3.3 (a) For a free particle of mass m in one dimension, Schrodinger’s equation takes the form 


pout) _. P Oat) 
Ot ——— 2m ðr? ` 
If Yı(x,t) and Y2(x, t) are both solutions of this form of Schrödinger’s equation, and cı and cz are any constants, 


we have 


(ing ges ) (c¥i (2,2) + ¢2W2(c,t)) 


ðt 2m Ox? 
z0 Tee OI OF 
= C1 (in + E) P(x, t) + C2 (ing + Li) Po(x, t) 


=c, X0+cx0=0, 


because Y; (x,t) and W2(x,t) separately satisfy the free-particle form of Schrédinger’s equation. So any linear 
combination of the form c, 1(x,t) + c2¥2(x, t) is also a solution of the free-particle form of Schrodinger’s 
equation, consistent with the principle of superposition. 
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(b) If the solution of Schrodinger’s equation corresponding to passage through one slit is Y4 and the solution 
corresponding to passage through the other slit is W32, the principle of superposition tells us that a normalized 
combination A(WV, + Y2) is also a solution of Schrédinger’s equation. The probability density for electrons 
arriving on the screen is given by |W, + W.|?, and this displays peaks where Y4 and Wy interfere constructively 
and troughs where Y; and Yə interfere destructively. 


Comments (1) An equation is linear if it can be written in the form L f(x), where Lis a linear operator. That is 
why we moved the right-hand side of Schrodinger’s equation to the left, and then used the method of Exercise 2.1. 


(2) The question refers to Schrodinger’s equation which should never be confused with the time-independent 
Schrodinger equation. Both equations depend on the system under study, but they have very different roles: 
Schrodinger’s equation tells us how any wave function for the system evolves in time, while the time-independent 
Schrodinger equation tells us about the possible energies of the system and the corresponding energy 
eigenfunctions. 


Topic 4 — Separation of variables 


Solution 4.1 (a) Substituting (x,t) = Y(x)T(t) in Schrodinger’s equation gives 


dT h? dy 
)ih— = —T(t)——> T(t). 
Wih E = -TOE + V@Hla)TW) 
Dividing through by %(x)T (t) and rearranging then gives 
2 g2 
: poe ae OY ray: 


T(t) dt p(x) 2m dz? 


The two sides of this equation depend on different variables, so they must both be equal to the same constant, E. 
We therefore get two equations: 


1 fe dy 
“Gamer SA 
1 . dT 
m u = E, 


which can be written as 


- 5S + V (2) (a) = Eyle) 


The general solution of the second equation is T(t) = Ce—iFt/h, where C is an arbitrary constant. (This constant 
is generally chosen to be equal to 1; there is no loss in generality in doing this because the arbitrary constants 
accompanying T(t) and y(x) are multiplied together in the product function ~(x)T(t).) 


(b) Not all solutions of Schrodinger’s equation can be written in the form of a product of a function of position 
and a function of time. Because Schrodinger’s equation is linear, any linear combination of solutions is also a 
solution, so a solution of the form 


W(z,t) = 5 ai pila) e/a 


a 


will satisfy Schrodinger’s equation, although it is not itself in a separated (i.e. stationary state) form. 


Comments (1) The natural way to start this question is to substitute the given form of wave function 

(x,t) = ¥(x)T(t) into Schrodinger’s equation. The key step is then to rearrange the equation so that one side 
involves only functions of x and the other side involves only functions of t. This is achieved by dividing by 
Y(2)T(t). 

(2) Note that the first line of the solution uses ordinary derivatives rather than partial derivatives. This is 


appropriate because y(x) and T(t) are both functions of a single variable, so 07/Ox and OT /Ot can only mean 
dy/dx and dT/dt, respectively. 
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Solution 4.2 Substituting y(x, y, z) = X(x)Y (y)Z(z) in the time-independent Schrödinger equation gives 
q’ X gY d?Z 

-E (YZ + xz + roro 
Dividing both sides by X(x)Y (y) Z(z) then gives 

can re ee ee ae) 1 PZ\_ 

2m \ X(x) dr? = Y(y) dy?  Z(z) dz? 
The three terms in the round brackets on the left-hand side are functions of different variables. The only way the 
left-hand side of the equation can be equal to the constant on the right-hand side is for each of these terms to be 
equal to a constant. We therefore have 

1 X _ 

X(x) dz? 7 


) = exer WZC). 


Z(z) dz 
Substituting these equations back into the unseparated equation, we see that 
h2 


Ag A = E. 
z— (Or + Cy + Cr) 


Comments (1) The method of separation of variables works for partial differential equations involving many 
variables, just as well as for those involving two variables. It can be applied step-by-step, splitting off one variable 
at a time, but here we have used it more efficiently relying on the fact that if a sum like f(x) + gly) + h(z) is 
equal to a constant for all values of x, y and z, it follows that f(a), g(y) and h(z) are each equal to constants. 


(2) In the last line of the solution, it turns out that Cy is negative, and is normally written as —k?, with similar 
notations for Cy and C’,. However, the proof of this fact requires consideration of the boundary conditions, and is 
not part of this question. 


Topic 5 — Born’s rule 
Solution 5.1 (a) In one dimension, Born’s rule states that the probability at time t of finding the particle in a 
small interval of width dx, centred on z, is |W (x, t)|? ôx. 
(b) The required probability is 
probability = |U (x, t)|? dx 


= (VF swam = (V2 sin(0.257) own) x 0.002L 


9 . . 
ave sin? (0.257) etiFt/he-iFt/h x 0.002L 


2 1 
=-= X= .002L = 0.002. 
T*7* 0.00 0.00 
This probability does not depend on time because the time-dependent phase factor en 


the modulus of the wave function. 


is removed by taking 


Comments (1) The wording of this exercise includes the words “small interval’. This is effectively an invitation 
to use Born’s rule without integrating over an extended range: we are assuming that the wave function 

is practically constant over the small interval from 0.249 to 0.2512, and that the appropriate value to use 
throughout this region is Y (0.25, t). 


(2) The probability works out to be a real number between 0 and 1, as it must. 
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Solution 5.2 The probability of finding the particle between x = 0 and x = L/5 at t = 0 is 


ee 5 
probability = T if sin? (=) dr: 


We change the variable of integration to u = 57a/L. Then x = 0 corresponds to u = 0 and x = L/5 corresponds 
to u = T. Also, x = (L/5a)u and dx = (L/5z) du, so 


2/L s 
probability = (=) f sin? udu 
ot J Jo 
L 


L 
2 
L 


5r 


Comments (1) In this exercise, the given interval cannot be treated as small in the context of using Born’s rule 
(no matter what the value of L). This is because the wave function varies significantly across the interval. We 
cannot use Born’s rule directly, but must integrate the probability density |Y (x, 0)|? over the interval, from x = 0 
to x = L/5. 


(2) This question could also be answered by drawing a sketch of |W(x,0)|? between x = 0 and x = L, noting 
that it consists of 5 peaks, each of the same area, and remembering that the area under the whole graph must be 
equal to 1. 
Solution 5.3 (a) For x > L/2, the probability density for finding the particle at position x is 

|ab1 (x) |? = (Ce-*)* Ce = Oj e7, 


Hence the probability of finding the particle anywhere with z > L/2 is 


eo 2ax T=00 —aL 
—2a | 


Poin = f Pe ar = of |E = 
A Life a=L/2 2a 


For x < —L/2, the probability density for finding the particle at position x is 
I(x)? = (Ce**)* Ce = |C Pe. 

Hence the probability of finding the particle anywhere with x < —L/2 is 

a S —aL 


—L/2 
Pecans | OPen jor E 


— 00 L=—0o 


The total probability of finding the particle somewhere outside the well is therefore 


Ic? 
Pout = Pestj2 + Pr<-L/2 = ae 


(b) The particle must be somewhere, so the probability of finding it to be somewhere inside the well is 


Cl? 
ee ee |f -oL 


Comments (1) The region outside the well consists of two subregions: —o0 < x < —L/2 and L/2 < x < œ. 
You need to include both of these to find the probability of detecting the particle outside the well. 


(2) Rather than evaluating the integral from —oo to — L/2 directly, you could change the variable in this integral 
from x to u = —; this would give an integral of the same form as that for the region from L/2 < x < ov, so the 
contributions from the two subregions are the same. 


(3) We have used the ground-state energy eigenfunction (rather than the corresponding stationary-state wave 
function) to find the probability because we know that the time-dependent phase factors in the stationary-state 
wave function Y and its complex conjugate W* cancel out in the product Y*W, leaving only a time-independent 
probability density (see Exercise 5.1 for example). The eigenfunction is equal to the stationary-state wave function 
att = 0. 
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Solution 5.4 The particle is most likely to be found at the point where |W(z,0)|? has its maximum value. We 
have 


|W (a, 0)|2 = |Al2a2e~227/0", 


The stationary values of this function are found by setting 


i= < (Ape) 


n |A 2 ($ x e720? /a? +r? x 2 (n) 
T T 


= |A 2 (2eer™"/ _ peri) 


a 


2 
= 2r Ar (1 _ — e722? /a? 
a 


This equation has solutions at x = 0 and z = +a/ v2. 


A rough sketch of the function shows that z = 0 corresponds to a minimum, and z = +a/./2 correspond to 
maxima, so the particle is most likely to be found at x = +a/ v2. 


Comment Some care is needed in doing the differentiation. The probability density involves a product of x? and 
e722°/a? which is differentiated using the product rule. The term e722°/a? must be differentiated using the chain 
rule: we differentiate the exponential function (this gives the same function back again) and also multiply by the 
derivative of its contents: 


ae tle! = eo 22*/a” x L (—20?/02) = 88 e, 


Solution 5.5 (a) In one dimension, Born’s rule states that the probability of finding the particle in a small 
interval of width ôx, centred on z, is |Y (x, t)|? da. Probabilities are dimensionless pure numbers, so |W (x, t) 


must have the dimensions of length’. Hence W(x, t) has the dimensions of length / ? and its SI units are 
-1/2 
m=, 


| 2 


b) 0 f(x) = eT" /2a” is unsuitable because it is dimensionless (it would have to be multiplied by a suitable 
factor with dimensions of length !/ ? in order to acquire the right dimensions. 

Gi) g(a) = (\/7/a) e~*’/2@ is appropriate on dimensional grounds because it has the dimensions of length~!/?. 
(iii) A(x) = (\/7/a) e~*”/? is unsuitable because the term e~*’/? has no consistent units. If it were expanded as 
a power series, it would give 


T =g? pa =z? P 
2 2! 2 Ea 


Since x is a position coordinate with the dimensions of length, successive terms have different units; such terms 
cannot sensibly be added together. To make physical sense, arguments of functions like the exponential function 
must always be dimensionless. 

Comments (1) It is worth remembering that any one-dimensional wave function has units of length !/ 2 and 
this can be a useful check. Note, however, that this is only true in one dimension. In three dimensions, wave 
functions have the dimensions of length—°/ 2; 


(2) Part (iii) relies on the convention used throughout physics where the position coordinate x is not a pure 
number, but is a number times a unit of length. This differs from some usages in maths, where the symbol x could 
perfectly well be a pure number (perhaps a ratio of two lengths). 
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Topic 6 — Normalization and phase 


Solution 6.1 For the given wave function, the normalization condition is 


boy )Y(0,0)P ax = [AP f et le de. 


oO 


Using the standard integral given in the question, this gives 


1 = |A| av7. 


A suitable choice for A is therefore 


ae () 


Any other value of the form 


1/4 
A= 1 e* 
Taz 


where a is real, could be chosen because any overall phase factor multiplying a wave function has no physical 
consequences. 


Comment In spite of the freedom noted above, we would normally set et% = 1, giving the real positive value of 
A originally calculated. This has the advantage of simplifying any subsequent calculations made with the wave 
function. 


Solution 6.2 We multiply the unnormalized eigenfunction (x) by a constant N, and then choose N to ensure 
normalization. Let 


v(x) = No(a). 


Then we require that 


1= f WPa 


—oco 


= |N}? f jġ(x)|? dx 


We change the variable of integration to u = 2rx/L. Then x = 0 corresponds to u = 0 and x = L corresponds to 
u = 2r. Also, x = (L/27)u and dx = (L/2z) du. So, using the standard integral given in the question, 


L 2m 
1 = |N? — sin? udu 


=|N|*— xa 


= |NŽ. 


Hence a possible choice of N is ,\/2/L and the corresponding normalized eigenfunction is 


(a= Žsn (=) . 


Comments (1) Given an unnormalized wave function, with no adjustable constant, you are free to multiply it 
by an (as yet) undetermined constant N. Then you can use the same procedure as in Exercise 6.1 to find an 
appropriate value for that constant. 


(2) Itis often convenient to start a normalization calculation with the number 1 on the left-hand side of the 
equation, allowing a chain of equalities to flow in subsequent steps. 


(3) Note how the limits of integration change from x = +00 to x = 0 and x = L because the integrand vanishes 
outside this region. 
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Solution 6.3 Normalization requires that 


=f ff Wey2)Pardyae 
L pl pL 2 
= lap f I f sin? (=) sin? (=) sin? (=) da dy dz 
D À L 2 i 
= ap f sin? (=) ae | sin? (=) a f sin? (=) dż. 


For the integral over x, we change the variable of integration to u = mg/L, for the integral over y, we change 
the variable of integration to v = 27y/L and for the integral over z, we change the variable of integration to 
w = 3nz/L. Then 


LLL T 2m 37 
1=]A? 2 — — sin? udu f sin? vdv f sin? w dw 
T 27 37 0 0 0 
o 2L L L nan 3r 
7 n 2m 3r 22 2 
L8 
-_ 2 
= |A aa 


so a suitable choice of A is (2/L)?/?. 


Comment Normalization in three dimensions involves a volume integral of the probability density. The present 
case uses Cartesian coordinates and, because of the product form of the wave function, the volume integral splits 
into a product of three separate integrals over x, y and z. 


Solution 6.4 we, We, Wp, and Yg all represent the same state as Ya because they are equal to Ya times a phase 
factor. We have: 


Ue = ida = tb 
be = —iva = e™ ypa 
bp = ea 

ig = T/a. 


Comment To spot whether ġa = a1 + agWo+... (with ay Æ 0) is a multiple of dg = bi, + boye+..., we 
can multiply ġa by the factor bı /a; that ensures its first term is the same as the first term in dg, and we can then 
check whether the remaining terms are also equal. 


Topic 7 — Square wells 


Solution 7.1 The energy levels for a particle of mass m in a three-dimensional cubic infinite square well with 
sides of length L are: 


nh? 
2mL? 
The ground state has quantum numbers ng = 1, ny = 1, nz = 1, and energy 
307h? 
2mL?- 
When the sides of the box are shrunk to L/2, the ground-state energy becomes 
Br? A? 12m? h? 
Qm(L/2)2  2mL ` 


(n2 + n3 + n2). 


Einitial = 


Efinal = 
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The change in energy is 
On hh? 
AE = Efna — Finiti = mE 


Comment The compression of a small quantum system can require considerable energy. This helps to explain 
why atoms display considerable rigidity in spite of being composed of very insubstantial particles (electrons) 
outside a tiny central nucleus. 


Solution 7.2 The energy levels for a particle of mass m in a three-dimensional cubic infinite square well with 
sides of length L are: 


2 


2mL2 


(nz + ny +72). 


The states that have energy 3h?/mL? have n2 + n? + n? = 6. The following combinations of quantum numbers 
achieve this: 


Nng =l, ny=1, nz=2, 
so the energy level is three-fold degenerate or its degree of degeneracy is 3. 


Comment Note that, in specifying a quantum state, it matters which quantum number has which value. 


Topic 8 — Continuity boundary conditions 


Solution 8.1 y(x) and dy/dx are continuous functions at any point where the potential energy function V (æ) is 
finite, including the walls of a one-dimensional finite square well. 


At the walls of a finite square well, y(x) and dy/dx are both continuous; at the walls of an infinite square well, 
w(x) is continuous but dy/dx is not. 


Comments (1) Asa finite square well becomes deeper and deeper, dy/dx changes more and more rapidly in the 
vicinity of the walls of the well. In the limiting case of an infinite well, dy/dx changes abruptly at the walls, 
corresponding to a discontinuity in the slope of the eigenfunction. 


(2) These boundary conditions apply to energy eigenfunctions; they do not apply to all wave functions (see p. 69 
of Book 1). For example, a wave function could have a discontinuous top-hat shape, as in Exercise 9.3 below. 


Solution 8.2 We are told that Yin = Acos(kx) and You = Ce are the forms of the eigenfunction inside and 
outside the well respectively. Then continuity boundary conditions require that 


dwvin out 
Vin(L/2) = Pour(L/2) an dz IL/2 dx |L/2 
The first condition gives 
A cos (=) = Ce ee, 
Differentiating the eigenfunction, 
dyi d l 
an = —kAsin(kx) and “ow = —aCe *, 


so the second condition gives 


—kAsin (=) = —aCe %#/?, 
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Dividing the equation expressing the second boundary condition by that expressing the first, and cancelling a 
factor of —1, we conclude that 


E 
k tan (=) =a. 


Comment As the square well becomes very deep, a becomes very large, and the eigenfunction decays very 
rapidly outside the well. In the limiting case of an infinitely deep well, œ tends to infinity and the above equation is 
satisfied by kL /2 = 1/2, that is by kL = 7, as expected. 


Solution 8.3 Denote the interval —L/2 < x < 0 by Region 1, the interval 0 < x < L/2 by Region 2 and 
elsewhere by Region 3. Then the time-independent Schrodinger equation in Region 1 is 


h? d?y(2) 
—— =F ; 
2m dr? Hiel 
Using the form of the given eigenfunction appropriate for Region 1, this gives 
R? q? 
-mT (Acos(kıx) + Bsin(kix)) = E(Acos(kıx) + B sin(kıx)). 
Carrying out the differentiations, we obtain 
hk? 
m (Acos(kix) + Bsin(kıx)) = E(Acos(kiz) + B sin(kız)), 
m 
so 
E= ake 
2m ` 


In Region 2, the time-independent Schrödinger equation takes the form 


h? dy 
-E THO) + You(e) = Eyle), 


and a similar argument gives 


Equating these two expressions for the energy eigenvalue then gives 


hike he kee 
A -aR yy 
2m 2m ay 


which is the required relationship between k1, k2, Vo and the mass m of the particle. 


Comments (1) Note that we can apply the time-independent Schrodinger equation separately in Region 1 
and Region 2, and our two expressions for the energy E can be equated because they refer to the same energy 
eigenvalue, E. 


(2) Ina sense, our final equation expresses basic energy accountancy: there is no potential energy in Region 1, so 
the kinetic energy in Region 1 (h?k?/2m) is equal to the sum of the kinetic and potential energies in Region 2 
(h?k3 /2m + Vo). However, this is talking rather loosely. The energy eigenvalue E refers to an energy 
eigenfunction which spreads throughout both Region 1 and Region 2. The particle is not confined to either 
Region 1 or Region 2; if it were, it would effectively be in a smaller well and its energy would be different! 


Solution 8.4 The eigenfunction must be continuous at x = —L/2 (the boundary between Regions 1 and 3 as 
defined in the preceding solution) so 


—kıL —kı L 
Acos ( > ) + Basin ( 5 )=0 


which simplifies to 


tan (BE) A 
an 5) = 
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The eigenfunction must also be continuous at x = L/2 (the boundary between Regions 2 and 3) so 
koL koL 
C cos (=) + Dsin (=) = (0. 


which simplifies to 


tan (PE) E 
an 5) = D 


At x = 0 (the boundary between Regions 1 and 2) both the eigenfunction and its derivative must be continuous, so 
Acos(0) + Bsin(0) = C cos(0) + D sin(0) 

and 
—kı Asin(0) + kı B cos(0) = —k2C sin(0) + k2 D cos(0) 

giving 
A=C and kB = kD. 


Combining the results obtained from all four continuity boundary conditions we conclude that 


an (EL) C RA kon (hk 
alao aB a ON)’ 


which is the required relationship between k1, k2 and L. 


Alternatively, we could write this as 


koL ky kL 
cot | —— ] = ——cot | —], 
2 k2 2 


Comment [In the limiting case of an infinite square well with no step in it, we have kı = k2 = k. From the first 
form of our final answer, we then obtain tan(kL/2) = 0, and from the second form we obtain cot(kL/2) = 0. 
These are equivalent to sin(kL/2) = 0 and cos(kL/2) = 0, which are the appropriate conditions for the odd and 
even energy eigenfunctions (respectively) of a symmetrical infinite square well. 


Topic 9 — The overlap rule 


Solution 9.1 (a) The probability p; of measuring the energy eigenvalue F; in a state described by the wave 
function (x,t) is 
2 


Pi = E (x)(x, t)dz| , 


where 7;() is the energy eigenfunction of the system corresponding to the eigenvalue Æ;. Both Y;(x) and U(z, t) 
are assumed to be normalized. 


(b) If y(x) is odd and W(x, t) is even, the integrand in the above integral is odd. This is integrated over 
a symmetric range, centred on the origin, so the integral is equal to zero. The probability of getting the 
corresponding eigenvalue Æ; is then equal to zero. 


Comment Note how the overlap rule works for energy measurements. We are looking for the probability of 
getting a given energy eigenvalue F; when the energy of a system is measured. We need to know: (i) the wave 
function U(x, t) that describes the state of the system at the time of measurement and (ii) the energy eigenfunction 
w;(a) that corresponds to the energy eigenvalue E;. This allows us to construct the overlap integral representing 
the probability amplitude for this measurement result. Don’t forget to take the square of the modulus of the 
probability amplitude to find the corresponding probability. 
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Solution 9.2 If the system is in a state described by the energy eigenfunction 7 (x), an energy measurement is 
certain to yield the corresponding eigenvalue, Æ. Consequently, the probability of obtaining any other eigenvalue, 
such as Fo, is equal to zero. 


According to the overlap rule, the probability of obtaining the eigenvalue E» in the state 71 (æ) is 
2 


OO 
EZ OLOLOE 
00 
which we know is equal to zero. Since the equation |z|? = 0 implies that z = 0, we conclude that 


[Slot ae = 0, 


confirming that yı (a) and w2(a) are orthogonal functions. 


probability = 


Comment  Self-consistency of the overlap rule requires energy eigenfunctions with different eigenvalues to be 
orthogonal. However, the energy eigenfunctions are solutions of the time-independent Schrodinger equation, 
Hu = Ew, which is the eigenvalue equation for the Hamiltonian operator, H. This implies that there must be 
something special about the Hamiltonian operator in order for it to deliver eigenfunctions with the required 
orthogonal property. In Book 2, you will see that this special something is called the Hermitian property of the 
Hamiltonian operator, named after the French mathematician Charles Hermite. 


Solution 9.3 The energy eigenfunction corresponding to the ground-state eigenvalue E1 = 1°h? /2mL? is 


ypı(z) = i Fe0s (=) ‘ 


so the probability of getting the value F; is 


lee) 2 
pi=|f via) ¥(x,0) dx 
i 
= te =) )yZ T 
w 
w/2 2 
= VIL F cos ( ) ax 
w/2 


Rather than change variables in the integral, we can directly use the fact that 


2 


f cos(kx) dz = 7 sin(ka) + arbitrary constant 


to obtain 
2 L l : S 
= |,/— — |sin (| — 
m wL 7 L/\—w/2 
2 
L TW 
= py Zan 
8L . 9 =) 
= —— siní | — 
T? w 2L 
Comments (1) When substituting in the limits of integration, we have used the fact that, if f(x) = —f(—2) is 


an odd function, then 
ro] = 27(0). 


(2) The ratio w/L is dimensionless, so the final answer is dimensionless, as required for a probability. 


(3) We can also check some values: at w = L, pı = 8/7?, and as w — 0, sin? (nw/2L) — (nw/2L)? so 
pı — 0. These values are reasonable because they both lie in the range between 0 and 1. 
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Topic 10 — The coefficient rule 


Solution 10.1 In order for the wave function to be normalized, we must have 


112 
1=|- k 
f + lal 
so 
3 
2n & 
|a] ~ 4? 


and a possible value of a is /3/2. 


This value is not unique and any choice of the form /3e'~/2, where a is real, would also produce a normalized 
wave function. In this case, however, different phases of a lead to different wave functions, representing different 
states. For example, 


V3 


1 3 1 
5 Vi(2) + 3 pala) represents a different state from 5Vi(2) — -z vale), 


with a different probability distribution for position measurements. 


Comments (1) In a situation like this, normalization can be achieved using the coefficient rule directly; there is 
no need to write down or evaluate explicit normalization integrals. 


(2) Distinguish this case (where the relative phase of two terms in a wave function matters) from the fact that an 
overall phase factor, multiplying the whole wave function, has no physical significance. 


Solution 10.2 Normalization requires that 
1 = |2N]? + |— V5N/? = 9N’, 
so we can choose N = 1/3. 


The probability of getting energy Æ; is 


Pi =z 


and the probability of getting energy Es is 


2 
=;/5 B 
zais 


5 


p3 = g 


The expectation value of the energy is 
5 14 


4 
(E) = p1 Ey + poE2 = 9° + 92€ = 9° 


Comments (1) We can check that pı + p2 = 5/9 + 4/9 = 1, in agreement with the normalization rule for 
probability. 


(2) In this case, the expectation value is not equal to either of the energy eigenvalues, but lies somewhere 
between them. 


Solution 10.3 The normalization condition is 

1 = | Ae 24/4)? + | — ideit]? — 9) Al? 
so a suitable value of A is 1/2. 
The probability of measuring energy F is 


— il —ik\t/h = ae 1 
Pl B f2 e = 2 , 
and the probability of measuring energy E» is 


—i | 
e 7 iBat/h 


p3 = V2 5? 
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so the expectation value of the energy is (EZ) = (Fy + E2)/2. 


Comment The fact that this expectation value is independent of time, even though the wave function and other 
expectation values depend on time, is related to the principle of conservation of energy. 


Solution 10.4 (a) The normalization condition is 


f 2 ; 2 F 2 
Jae irh + Jaze "| + aget] — lal? + |a9|? ae \a3|? = 4 


(b) No. It never makes a difference to multiply a whole wave function by a phase factor (a complex number of 
unit modulus). 
(c) Yes. This will affect the interference between different terms in Y (x, t). 
(d) No. The expectation value of the energy is 
(E) = |a1|?E1 + |az|? £2 + las)? Es, 
which is unchanged by reversing the sign of a1. 
Comment Reversing the sign of a; alone changes the state of the system (witnessed, for example, by the change 


in the expectation value of x) but the original state and the new state share the probability distribution for energy. 
So, in general, the probability distribution for energy is insufficient to specify a quantum state. 


Topic 11 — The sandwich integral rule 


Solution 11.1 (a) The sandwich integral rule states that, for a one-dimensional system described by the 
normalized wave function (x,t), the expectation value at time t of an observable A is 


(A) = a W* (x,t) AV(a, t) dz, 


where A is the quantum-mechanical operator corresponding to the observable A. 
(b) The expectation value of position is 


(a) =| Ate-°” Ae? dx 


=00. 


OO 
= a f ve 2" dg 


o0 
= 0, 

because the integrand is an odd function and the range of integration is symmetrical about x = 0. 

The expectation value of momentum is 


(Px) =] Loe (-ing aero) da 
g 


—oo 


= —ih|A|? L eon” (—2are-**" ) dx 


—Cco 
2a 2 
= 2inal AP f ce Oe dg 
=00 
=0, 


because the integrand is odd and the integral is over a range centred on x = 0. 


Comments (1) When using the sandwich integral rule, allow the relevant operator to act on the wave function 
Y (x,t) before multiplying by the complex conjugate Y* (x, t), collecting terms and integrating. This does not 
matter much in part (a) but is more important in part (b). 


(2) The integral in part (b) can be evaluated more simply by noting that —ihO/Oz is an odd operator (changing 


sign when x — —2) while e~°” is an even function. This shows that the integrand is an odd function (without the 
need to carry out any differentiation). 
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Solution 11.2 We have 
(0?) = | WE (a) iii 
9 L 


= al x sin? (=) dz. 


Change variables to u = nra/L. Then x = 0 corresponds to u = 0 and x = L corresponds to u = nz. Also, 
x = (L/nz) u and dz = (L/nz) du, so 


27 Le 
(2?) = = (=) f u? sin? udu. 
L \ nr 0 


Using the standard integral given in the question, we conclude that 


2(L\e% nn 
2a Ff |) M8 (9,272 _ 
(x == (=) 19 (2” n“ — 3) 


I? 3 
= — | 2— ; 
6 ner? 
Comments (1) The answer is dimensionally correct because the expectation value of x? must have the units of 
length”. 


(2) In the limit of very high quantum numbers (n — 00) we get (x?) — L?/3. Compare this with the 
corresponding classical result. In classical physics, the particle is equally likely to be found anywhere inside the 
box, so the probability density has the constant value 1/L for 0 < x < L, giving a classical expectation value for 
x? of i. x?(1/L) dx = L?/3. In the limit of very high quantum numbers, the quantum result approaches the 
classical result; this agrees with the correspondence principle (Book 1 p. 135). 


Solution 11.3 In the given state, the expectation value of p? is 
OO 
(2) = =s fete pte Pe ae 
—oco 


TA 


h? i — x? /2a? a —2? /2a? 
=e, fee Ta (7) ae. 


Using the derivative given in the question, we obtain 


2 oo 2 
o=- fe (G1) ee ae 
Tta J a? \a 


We change variables to u = x/a and use the standard integrals given in the question. Then z = au and dx = a du. 
The limits of integration remain oo, so 


2 oo lo) 
(po) = -57 ef uet dua f ew au] 
he [yT 
ae $ - v7 
h2 
~ 2a?’ 


Comment As the constant a becomes smaller, the wave function becomes more narrowly peaked around 

x = 0, and our answer shows that (p2) increases. This can be understood as follows. When the wave function is 
represented as a linear combination of plane waves, a greater spread of k values is needed to produce a narrow 
peak rather than a broad peak. This corresponds to a greater spread in momentum values py = ñik, and so to a 
greater expectation value of p2. The greater spread of momentum values is also consistent with the Heisenberg 
uncertainty principle. 
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Solution 11.4 The expectation value is given by 


2,2 1 CO 2 )52 _ 4222 nion 
(e a ms x“ /2a e ME e x*/2a är. 
vra J- 


We can collect together all the exponential factors to obtain 


=A o l i — (0241/02)? 
(e )= Waa [o dz. 


Using the standard integral given in the question, with 


Tlia? 
a=) +1/@ = UPAR z2 2 
we get 
An 1 Ta _ 
(e : ) yra * (1 te (1 Na") ve. 


Comments (1) The expectation value of any function of x could be calculated in a similar way. 


(2) One way of checking this answer is to consider the limiting case A = 0. In this case, en” = 1, so the 
sandwich integral becomes equal to the normalization integral, which is equal to 1. Our final answer gives this 
value for \ = 0. 


Topic 12 — Uncertainties and the uncertainty principle 
Solution 12.1 Energy E; = € occurs with probability 


Pi =l 


Energy E> = 2e occurs with probability 


Jij? 1 
P2 = 3 FF 
The expectation value of the energy is therefore 
3 1 5 
E) = — -2e = —€. 
(E) ze t qe = Ze 
The expectation value of the square of the energy is 
3 1 7 
Py a PP foe a 
(E?) = Že? + 7 (2e)? = $e 


So 


(AE)? = (E?) - (E) 
‘2 25 9 2 


4 16 I6 
and the uncertainty of the energy is 
V3 


AE = —e. 
a 


Comments (1) The probabilities are found by taking the square of the modulus of the coefficients in the wave 
function. These probabilities turn out to lie in the range from 0 to 1, as required. 


(2) The expectation value of E is a weighted sum of the possible values of Æ, with each allowed energy value 
multiplied by its corresponding probability. The expectation value of Æ? can be written as a similar sum with the 
same probabilities, which now multiply the squares of the values of E. This works because, to find the square of 
the energy in any particular measurement, we can measure the energy value E and then take its square. 


(3) It is convenient to calculate the square of the uncertainty first, and then take the square root. 
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Solution 12.2 (a) The square of the uncertainty in « is 
(Ax)? = (x?) — (x)? = 2.0 nm? — (1.0nm)? = 1.0nm”, 
so 


Az = 1.0nm. 


(b) The minimum value of Apg, consistent with the above value of Az is 
ñ _ 1.06 x 10°-**Js 
2Ax  2x10-9m 


The minimum Ap, consistent with the given value of Az is zero. 


(Aps)min = = 5.3 x10 ke ms. 


Comment The Heisenberg uncertainty principle applies to position and momentum components along the same 
axis. There is no uncertainty principle for a position component along one axis and a momentum component along 
a perpendicular axis. 


Topic 13 — Harmonic oscillators 


Solution 13.1 We have 


(e=) = 4 (—2an0-**") 


= (—2a + 4o2x?\e~0?” 


Substituting into the time-independent Schrodinger equation, we obtain 
h? 2,2) .—axr? 1 2,—a07 —aa? 
-z 2e + de ae + 3Ca*e% = Ee : 


2 . 
o”" and rearranging, we conclude that 


Aha? 1 hi 
(- . +50) e+ (#2 - £) =0. 
2m 2 m 


For this equation to be true for all x, each of the terms in round brackets must vanish separately, so we have 


2.2 2 
aut = lo and ri sF: 
2m 2 m 


Dividing through by e7 


Hence 


Comment Notice that two pieces of information can be extracted from a single equation. This is because the 
time-independent Schrodinger equation is valid for all values of x. Putting x = 0 gives a relationship between Æ 
and a. Using this relationship and considering the time-independent Schrödinger equation for x Æ 0 then gives a 
relationship between a and C. 


Solution 13.2 The energy levels of a simple harmonic oscillator are given by 
E= (n+ 5) hwo forn = 0,1,2,..., 


where wo is the classical angular frequency. The zero-point energy is the lowest possible energy above the bottom 
of the potential well: 


1 
Fzero-point = fo = zwo. 


The radiative transitions in a harmonic oscillator are subject to a selection rule that restricts An to +1. So a 
radiative transition from a state with n = 5 to a state of lower energy can only be to the state n = 4, and the energy 
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of the photon emitted is 
Ephoton = (5 F Z )ħwo = (4 F 5) hwo 
= hwo 


= 2 Ezero-point 
= 3.0eV. 
For electromagnetic radiation, Æ = hf and c = fA, so the wavelength of the emitted light is 
c he 663% 107°4J s x 3.00 x 108m s7! 


zts La LI O SS A i a 
F E 3.0 x 1.60 x 10-9 J A 


Comment Because of the selection rule, and the fact that the energy levels of a harmonic oscillator are equally 
spaced, a harmonic oscillator produces only one spectral line no matter how much it is excited. Harmonic 
oscillators are unique in this respect; atoms have much more complicated spectra. 


Solution 13.3 In terms of raising and lowering operators 


= = (A+4) and B= a hâ- Al). 


Combining these results and multiplying out the brackets (whilst preserving the order of operators) gives 


RP, = (A +4") (A- ât) 


ih 
— (AA — A'A! + ATA — AA") 
= = (AA - AAt 1) , 
where we have used the commutation relation AA! — ATA = 1. 
So, in a harmonic oscillator energy eigenstate with quantum number n, 


ep.) = f yle) (AA- ART — 1) yala) de. 


—oo 


This splits into three integrals. We have 
f @ RR eno) dr =a f yila) bala) dr =0, 
—oo —oo 


where œ is a proportionality constant. The integral is equal to 0 because n(x) and Yn+2(x) are orthogonal. 
For n = 0 or n = 1, we have AAdn (2) = 0, while for n > 2, 
Co OO 
f E Bun (a) ax = 8 [v5 (2) bn-ale) de = 0, 
=00 =00 
where (3 is another proportionality constant. Again, the integral is equal to zero because Yp (x) and Yn—2(x) are 


orthogonal. 


We therefore conclude that 


eb Ut SF g 
GB.) = +5 f vila) Yala) de 


“= 
because n(x) is normalized. 


Comments (1) When multiplying out brackets involving raising and lowering operators it is essential to preserve 
the order of the operators because they do not commute. 


(2) When arguing that a given integral involving raising and lowering operators is equal to zero, you need not 
give explicit values for proportionality factors such as a and (3 above. 


Solution 13.4 (a) Multiplying out the brackets (whilst preserving operator order in products) we obtain 
CEO re ee ere ee eee 
a ee oe aT 
+(AA-AA -AA+AA) 
Farner ed: 
=2(AA+A A). 


So the integral given in the question becomes 
j=% f wt (x) AAdn (x) da +2 f wt(x) A Alaa) de. 


Both terms on the right-hand side are equal to zero (by the argument given in the answer to the previous question). 
(b) Using the formulae given in the question, we find that, in the state Yn, 


+(e) = ; f 7 ya) (Â + A')an(x) de 


a2 


a2 


o= -3 | WoA- 3) unla) ae 


Subtracting these two equations then gives 


az oo D = a a 
ae) Fee) =5 | AA+ + A- 2h] yale) =0, 


a 


by the result derived in part (a). Hence (x?) /a? = (a?/h)(p2), as required. 
Using a* = h?/mC, we can rewrite this result as 


1 2 1,9 
5 C(t") = 5 Peh 


which shows that the expectation values of the potential and kinetic energies are equal in any energy eigenstate of 
a harmonic oscillator (i.e any stationary state of the oscillator, represented by an energy eigenfunction). 


Comment This quantum-mechanical result is closely related to a similar result in classical physics. For a 
classical harmonic oscillator, the time-average of the kinetic energy is equal to the time-average of the potential 
energy when both are averaged over a whole number of cycles. 


Topic 14 — Bound-state wave packets 


Solution 14.1 (a) At time ¢, the wave packet is 
1 ; ; 
U(z,t) = — ( r e`iEit/h i T pa] 
= a P o(2) 
-A 
E 


e- 3iwot/2 


= Ea (p(x) +iyp(x) a) ; 


(b) To see if there is a time when Ẹ (x,t) is orthogonal to the initial wave function, Y (x, 0), we consider 


(v(x) e Siwot/2 + i yox) a) 


oo e™3iwot/2 foo i ! 
/ ARE a a a f [va (x) + i ve(@)]" [ya (2) + i pole) e™o] da 


—oo —oo 


esiwot/2 eS | 
a a f [at (a) —ib3(x)] [di(w) + i pole) e™] da. 
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Multiplying out the brackets and remembering that energy eigenfunctions with different energy eigenvalues are 
orthogonal gives 


fore) e~ 3iwot/2 
[0,400 dr =F | f vile) ede + v5le) dala) eM] a 


—0o 
—3iwot/2 
= —— [1+e 0] , 
where we have used the fact that y (a) and w2(a) are normalized. 
The overlap integral between W(x, t) and (a, 0) is equal to zero when e~ “0! = —1. Since e”" = —1 for 
n = 1,3,5,..., the positive values of t for which the overlap integral is vanishes are: t = 1/wo, 37/wo,.... 


Comments (1) When an initial wave function U(x, 0) is given as a sum of energy eigenfunctions, the 
corresponding time-dependent wave function (x,t) is constructed by inserting appropriate time-dependent phase 
factors. Each energy eigenfunction has a different phase factor. Because these vary at different rates, the shape of 
the wave function can vary in an interesting way. 


(2) To establish orthogonality, we must evaluate an integral: the calculation in part (b) would not be valid 
without integral signs. 


Solution 14.2 (a) The probability density is 
1 ' ; . : 
W* (a, t) U(a,t) = 5 (h(a) AP + pala) 7n) (yola) ei + pafa) etn) 


= 5 (Ivot) + |abo(a) |? + (2) palæ) ei(Eo—Ea)t/h 


+ U5(2) volg) Eon), 


This is periodic with period T where (E> — Eo)T/h = 27. So 
= wh wh č r 
o E — Eo i 2hw E wo 


(b) wWo(a) and w(x) are harmonic oscillator energy eigenfunctions with n = 0 and n = 2; these are both even 
functions of x, so U(—x,t) = (x,t) at all times. Consequently, the probability density is an even function of x 
and the particle is always equally likely to be found in the regions x < 0 and x > 0. The motion of this wave 
packet does not correspond to a classical to-and-fro motion, but to a symmetric breathing in and out (in which the 
probability density undergoes one cycle of oscillation in half the classical period of oscillation). 


Comments (1) By contrast with the previous exercise, calculations of the probability density involve no 
integration. 


(2) The time-dependence of the probability density arises from the cross-product terms when we multiply out 
the brackets. For some wave packets, such as that considered here, the period of oscillation can be shorter than the 
classical period. 


Topic 15 — Free-particle wave packets 
Solution 15.1 At iis later time, we have 


U (a, t) t) i(ka—E;,t/h) dk, 


ay, 27 Fal 
where Ep = h?k?/2m. 
A(k) is the momentum amplitude function. The probability of finding the momentum to be in a small range of 


width fi dk, centred on hk is |A(k)|? ôk. More generally, the probability of finding the momentum to be in the 
range from hk, to fika is 


k2 
probability aj |A(k)|? dk. 
kı 
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Comments (1) The method used to find the time-dependent wave function is similar to that described in 
Exercise 14.1 for a harmonic-oscillator wave packet. The only difference is that we now have an integral over k 
rather than a discrete sum (and, of course, we must use the energy eigenvalues for a free particle, Ep = h2k?/2m). 


(2) The interpretation of A(k) is based on Born’s rule for momentum, which is closely analogous to Born’s rule 
(for position). Note, however, that A(x) is expressed as a function of the wave number k, while the momentum is 
Px = hk. This means you must remember to include various He of A in the descriptions of momenta and 


momentum ranges, but not in the formulae | A(k)|? ôk or Si Mes |A(k)|? dk. 


Solution 15.2 (a) The momentum amplitude function is 


1 oe : 
A(k) = Ta P(x, 0) e** da 
1 L/2 i f 
= a eiae e ike dz 
V2T J-L/2 VL 
L/2 


o 


-L/2 
—i(a—k)L/2 


where we have used the identity sin 0 = (e? — e~!9) /2i. 


(b) The probability of measuring a momentum between —fa and +ħa is 


+a 
probability = f |A(k)|? dk 
—a 
_ 2 L sin?[(a — k)L/2] dk 
tL (a — k)? l 
We change the variable of integration from k to u = (a — k)L/2. Then (a — k) = 2u/L and dk = —(2/L)du 
Also, k = —a corresponds to u = aL and k = +a corresponds to u = 0, so we have 


2 0 sin? u 2 
bability = — ——~ x —-d 
probability o 74” 


1 (@ sin? u 
= — 5 du. 
T 0 u 


For the special case a = 1/L, the standard integral in the question gives 


+. 2 
A gi es a GORE. 
T 


1 fl 
probability = — f 
T JO 


u2 


Comment The given wave function can be described as a truncated plane wave because it looks like a plane 
wave in the region between — L /2 and L/2 and vanishes outside this region. In spite of this fact, it corresponds to a 
spread of momenta ñk with a probability distribution given by |A(k)|?, which is nothing like the single momentum 
value associated with a plane wave that extends indefinitely. The reason for this is as follows. If the wave function 
is represented as a linear combination of e'** terms, a broad range of different k values is needed to get the 
required complete cancellation in the region |x| > L/2. The different k values correspond to different momenta. 
In accordance with the Heisenberg uncertainty principle, the momentum uncertainty increases as L decreases. 
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Solution 15.3 (a) The function Vyg(x, t) cannot be normalized. It corresponds to a particle with a definite 
value of momentum and an infinite uncertainty in position. Such a particle could be found with any value of x 
between —oo and +00, which is not realistic. 


(b) The probability that a momentum measurement will yield a value between pı = hk, and po = ko is 


k2 
/ | A(k) |? dk. 
ky 


The condition for normalization is that 


T | A(k)|? dk = 1. 


Comment In physical terms, the above normalization condition expresses the fact that a measurement of the 
momentum is bound to give some value between —oo and +00. A mathematical result known as Plancherel’s 
theorem (or Parseval’s theorem) guarantees that this condition also gives 


f Heta =i 


—co 


so the wave function is normalized. 


Topic 16 — Ehrenfest’s theorem 


Solution 16.1 Ehrenfest’s theorem states that 
d(x) _ (Px) 


dt m 


d(pz) La: ov 
dt Oz f° 
If V(x) = mgz, then OV/Ox = mg and 


d (px) 
dt 


and 


=—mg. 


Hence, 


d*(x) _ d (2) d (22) _ a mg 


d dt\ dt d m 7 


dt 


m 


Comment This is the quantum-mechanical version of the familiar classical statement that close to the Earth, in 

the absence of air-resistance, particles move with a constant downward acceleration of magnitude g. The quantum 
result differs from the classical one because it involves an expectation value. For a particle dropped from rest, the 

time of fall would have a probability distribution associated with the spatial spread of the particle’s wave packet. 


Topic 17 — Scattering and tunnelling 


Solution 17.1 Continuity of y(x) at x = 0 gives 
A+B=C. 

Continuity of dy/dx at x = 0 gives 
ikı(A — B) = ik2C. 


So, eliminating C, 
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Hence 
B(kı/kə + 1) = A(kı/k2 = 1) 


and 
B- kı — k2 


A ktk? 


Comment As a quick check, note that B/A = 0 when kı = kə. Sensibly enough, this means there is no 
reflection when there is no step. 
Solution 17.2 Continuity of y(x) at x = 0 gives 
A+B=C. 
Continuity of dy/dax at x = 0 gives 
ikı(A = B) = -ac. 
So, eliminating C, 
ik 
At p= AD, 
a 
Hence 


B(ikı/a — 1) = A(iki/a+ 1) 


and 
B ikı/a+1 ikı+a 
A ikja—l ikj—a 
So 
f ? (iki +a) . (iki + a) _ kta 
A| (iki a) ~ (~ika) kh? +02 


Comment |B/A)? is the reflection coefficient, so the physical interpretation of this result is that all the particles 
in the beam are eventually reflected. Particles can enter the region beyond the step (by barrier penetration), but 
they do not stay there forever and are eventually reflected. 


Solution 17.3 The probability current is 
ih „OY oy* 
s= -E (we we), 


2m Ox Ox 
We have 
y(x) = Ae? and y” (z) = Le 
a. =iAke*® and ii = —iA*ke™™*? 
Ox Ox 
so 
ja = "n (Ate (Ak elk) = Aci (—i A"k e=”) ) 


faeces 
= -ZA (ik + ik) 
= PE ap. 
m 


Comment In classical terms, the answer is the probability density of the beam, multiplied by the velocity of the 
beam, which is appropriate for a probability current. 
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